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We review the application of non-Abelian discrete groups to Tri-Bimaximal (TB) neutrino mixing, which is 
supported by experiment as a possible good first approximation to the data. After summarizing the motiva- 
tion and the formalism, we discuss specific models, mainly those based on A4 but also on other finite groups, 
and their phenomenological implications, including the extension to quarks. The recent measurements of 
#13 favour versions of these models where a suitable mechanism leads to corrections to #13 that can natu- 
rally be larger than those to #12 and #23- The virtues and the problems of TB mixing models are discussed, 
also in connection with lepton flavour violating processes, and the different approaches are compared. 
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1 Introduction 

Neutrino mixing [ 1-6 1 is important because it could in principle provide new clues for the understanding of 
the flavour problem. Even more so since neutrino mixing angles show a pattern that is completely different 
than that of quark mixing. The bulk of the data on neutrino oscillations are well described in terms of three 
active neutrinos. By now all three mixing angles have been measured, although with different levels of 
accuracy (see Tab. |T| ||7][8j). In particular, we have experimental evidence for a non vanishing value of the 
smallest angle #13 (see Tab. [2]|9 - 12|): considering the most precise results from DOUBLE CHOOZ, Daya 
Bay and RENO, we get 

sin 2 13 = 0.0253 ± 0.0035 , (1) 

for both the mass orderings. 

Models of neutrino mixing based on discrete flavour groups have received a lot of attention in recent 
years fl4f{T9) . There are a number of special mixing patterns that have been studied in that context. Most of 
these mixing matrices have sin 2 823 = 1/2, sin 2 6*13 = 0, values that are a good approximation to the data, 
and differ by the value of the solar angle sin 2 6*12. The observed sin 2 #12, the best measured mixing angle, 
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Table 1 Recent fits to neutrino oscillation data from (7|8). In the brackets the IH case. )* In this case the 
full (0, 2tt) is allowed. 



Quantity 



Fogli etal. [7| 



Schwetz et al. 181 

LL 



(10 



Am Ln 

AmL (10 



" 5 eV 2 ) 
" 3 eV 2 ) 



sin 2 6»i2 
sin 2 6> 23 
sin 2 6»i3 
Scp/tt 



' - 04 -0.22 

2.43±2-SS (2.42+ 07 ^ 



J -0.09 

S07+ 018 



-o.i 



3Q8+ 003 

u - oyo -0.026 



(0408l ;°l 5 ) 



0.0245 
0.89 



0.0034 m 0246 +u 
0.0031 tu.uzio_ 



+0. 0034s. 
.0031^ 



0.29 m qf) +0.32s 

-o.44 (u.yu_ , 43 ) 



7.62 ±0.19 



o ro+0.08 
z - do -0.10 



(2.40i ; ° 7 ) 



320+ 015 
0.491°-°* (0.53+ 05 ^ 



-0.05 

0.026+0-003 

1+0.54 



-0.07) 
0.004 <> u - uz ' -0.004- 1 



0-83±g;|| (0.07*) 



Table 2 The reactor angle measurements from the recent experiments T2K [9], MINOS 1 10 1, DOUBLE 
CHOOZ (11), Daya Bay (12) and RENO (B), for the normal (inverse) hierarchy. 



Quantity 



sin 2 26» 



13 



sin 2 0i3 



0.11+8:8 (0.14±g;i|) 0.028ig^J (0.038±85g) 



T2K [9 
MINOS 



DC [11 
DYB 
RENOTO 



Sj0| 0.041±8;^I (0.079l2:8g) 



o.oioi83S (o.Q2oi8:8i5) 



0.086 ± 0.041 ± 0.030 
0.092 ±0.016 ±0.005 
0.113 ±0.013 ±0.019 



022+ 019 
0.024 ±0.005 
0.029 ±0.006 



is very close, from below, to the so called Tri-Bimaximal (TB) value |20 -24) which is sin 2 0i2 = 1/3 (see 
Fig. |T|l. Alternatively it is also very close, from above, to the Golden Ratio (GR) value 1 25-28 1 which is 

0.276, where <f> = (1 + v5)/2 is the GR (for a different connection to the 



sin 



2 

5+V5 



GR, see Refs. |29 30)). On a different perspective, one has considered models with Bi-Maximal (BM) 
mixing, with siir ?i2 = 1/2, i.e. also maximal, as the neutrino mixing matrix before diagonalization of 
charged leptons. This is in line with the well-known empirical observation that 0i2 + 0c ~ 7r/4, where 
Be is the Cabibbo angle, a relation known as quark-lepton complementarity pl|j5T[ . Probably the exact 
complementarity relation becomes more plausible if replaced with 0i2 + O(0c) ~ t/4 (which we call 
"weak" complementarity). One can think of models where, because of a suitable symmetry, BM mixing 
holds in the neutrino sector at leading order and the necessary, rather large, corrective terms for 12 arise 
from the diagonalization of the charged lepton mass matrices ] 3 1 -57). These coincidences cannot all be 
relevant and perhaps all of them are pure accidents. But if one or the other of these coincidences is taken 
seriously then one is led to consider models where TB or GR or BM mixing are naturally predicted as a 
good first approximation. 

In the following we will concentrate on TB mixing which is perhaps the most plausible and certainly 
the most studied first approximation to the data. The simplest symmetry that, in leading order (LO), leads 
to TB is A4, the group of even permutations of 4 objects, a subgroup of S4 that includes all 4-object 
permutations. Thus, in the following, we will devote a special attention to A4 models, but alternative 
theories of TB mixing will also be briefly considered. The plan of the paper is as follows. In Sect. [2] we 
recall the definitions of TB, GR and BM mixing and the symmetries of the corresponding mass matrices. 
In Sect. [5] we summarize the group theory of A4. In Sect. [4] we review the structure of A4 models of lepton 
masses and mixings and, in two separate subsections, we first describe the baseline models and then those 
special models ] 58][59) where additional dynamical ingredients allow that the angle 0i3 can naturally be of 
different (and larger) order of magnitude than the deviations of 0i2 from the TB value. We also discuss the 
comparison with present data of the two options. In Sect. [5] we discuss the possible extension of the TB 
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i 1 — i 1 

sin" l2 5 + Vs i 2 

GR TB BM 

Fig. 1 The values of sin 2 612 for TB or GR or BM mixing are compared with the data at la. 

models to include quarks, possibly also in a GUT context. Our speculations on the origin of A 4 either as a 
subgroup of the modular group or as a remnant of an extra dimensional spacetime symmetry are presented 
in Sect. [6] A number of alternative theories of TB mixing are briefly considered in Sect. [7] Sect.|8]contains 
a summary on the implications for lepton flavour violation of the different models described in Ref. [ 58 1. 
Finally in Sect. [9] we derive our conclusions. 

2 Special Patterns of Neutrino Mixing 

Starting from the PNMS mixing matrix U (we refer the reader to Ref. (T|[5j for its general definition and 
parametrisation), the general form of the neutrino mass matrix, in terms of the (comple:^ mass eigenvalues 
m 1 ,m 2 ,m 3 , in the basis where charged leptons are diagonal, is given by 

m v = Z7*diag(mi, 1712, m%)w . (2) 

We present here some particularly relevant forms of U and m v that are important in the following. We start 
by the most general mass matrix that corresponds to 6\% = and 823 maximal, that is to U given by (in a 
particular phase convention) 

(C12 S12 \ 

~s 12 /V2 c 12 /V2 -l/y/2 , (3) 
S12/V2 C 12 /V2 l/y/2 J 

with C12 = cos 9 12 and s\ 2 = sin 6*12- By applying eq. |2} we obtain a matrix of the form 1 60 - 75) : 




(4) 



with complex coefficients x, y, z and w. This matrix is the most general one that is symmetric under 2-3 
(or fi — t) exchange or 

■m v = A 2i m v A 2Z , (5) 
where A 2 ^ is given by 




A 23 = 1 . (6) 



We absorb the Majorana phases in the mass eigenvalues m;, rather than in the mixing matrix U. The dependence on these 
phases drops in neutrino oscillations. 
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The solar mixing angle 6*12 is given by 

s - 2 2(9 = 8\x*y + y*(w + z)\* 8y 2 

12 8\x*y + y*(w + z)\ 2 + (\w + z\ 2 - |x| 2 ) 2 (a - w - z) 2 + 8y 2 ' W 

where the second equality applies to real parameters. Since #13 = 0, in this limit there is not no CP 
violation in neutrino oscillations, and the only physical phases are the Majorana ones, accounted for by the 
general case of complex parameters. We restrict here our consideration to real parameters. There are four 
of them in eq. Q which correspond to the three mass eigenvalues and one remaining mixing angle, 6*12. 
Models with fi-r symmetry have been extensively studied 1 60 - 77j . 



The particularly important case of TB mixing is obtained when sin 26»i2 = 8/9 or x + y = w + zf\ In 
this case the matrix m„ takes the form 

(8) 




In fact, in this case, U = Utb is given by [20-24) 



v/2/3 1/V3 

Utb = ( -1/V6 l/VS -l/y/2 ) . (9) 
-1/V6 1/V3 1/V2 

Note that Utb is a rotation matrix with special, fixed angles: indeed all the entries of Utb are pure 
numbers. This property is related to the particular pattern of m v which belongs to the category of the form 
diagonalizable mass matrices (78). These matrices are diagonalized by unitary transformations that are 
independent from the eigenvalues. At the LO discrete flavour models give rise to form diagonalizable mass 
matrices and the physical mixing angles are thus unrelated to masses. From eq. one obtains 

m^m! $1 $[ + m 2 $ 2 ^ + m3$3*J, (10) 

where 

$1 = -^(2, -1,-1), ^ = -^(1,1,1), -^(0,-1,1) (11) 

are the respective columns of Utb and m t are the neutrino mass eigenvalues. It is easy to see that the TB 



mass matrix in eqs. ( 10 1 and ( 1 1 1 is indeed of the form in eq. §8j. All patterns for the neutrino spectrum 
are in principle possible. For a hierarchical spectrum 7713 >> m-2 >> mi, to 2 ~ Am 2 tm , m 2 /m§ ~ 
^ m soi/ '^■ m atm an( l m i could be negligible. But also degenerate masses and inverse hierarchy can be 
reproduced: for example, by taking m 3 = — m 2 = n%i we have a degenerate model, while for m\ — — m 2 
and m 3 = an inverse hierarchy case is realized (stability under renormalization group running strongly 
prefers opposite signs for the first and the second eigenvalue which are related to solar oscillations and 
have the smallest mass squared splitting [42,79 - 84)). 

Note that the mass matrix for TB mixing, in the basis where charged leptons are diagonal, as given in 
eq. (|8j, can be specified as the most general matrix which is invariant under // — r (or 2-3) symmetry and, 
in addition, under the action of a unitary symmetric matrix Stb (actually S TB = 1 and [Stb, ^23] = 0): 

m y = Stbiti^Stb , m v = A 2 3m u A 2 3 , (12) 

where Stb is given by 

Stb =H 2 -1 2 | . (13) 




= w + 2 gives rise to TB mixing in another phase convention and is physically equivalent to 
: w + z. 
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Similarly, it is useful to consider the product m 2 = m\m e , where m e is the charged lepton mass matrix 
(defined as ip R m e ipi,), because this product transforms as m' 2 = UJm 2 U e , with U e the unitary matrix 
that rotates the left-handed (LH) charged lepton fields. The most general diagonal m 2 is invariant under a 
diagonal phase matrix with 3 different phase factors, 

m\ m e = T' rn\ m e T , (14) 

and conversely a matrix m|m e satisfying the above requirement is diagonal. If T n = 1 the matrix T 
generates a cyclic group Z n . In the simplest case n = 3 and we get Zj, but n > 3 is equally possible. In 
the n = 3 case we have 

/ 1 \ 

T TB = u . (15) 
V u? j 

where uj 3 = 1, so that T TB = 1. 

We are now in a position to explain the role of finite groups and to formulate the general strategy to 
obtain the special mass matrix of TB mixing. We must find a group G / which, for simplicity, must be as 
small as possible but large enough to contain the S and T transformations. A limited number of products 
of S and T close a finite group G t . Hence the group G / contains the subgroups Gs and Gt generated 
by monomials in S and T, respectively. We assume that the theory is invariant under the spontaneously 
broken symmetry described by Gf. Then we must arrange a breaking of Gt such that, at leading order, 
Gf is broken down to Gs in the neutrino mass sector and down to Gt in the charged lepton mass sector. 
In a good model this step must be realized in a natural way as a consequence of the stated basic principles, 
and not put in by hand. The symmetry under A23 in some cases is also part of Gf (this the case of 54, the 
permutation group of 4 objects) and then must be preserved in the neutrino sector along with S by the Gf 
breaking or it could arise as a consequence of a special feature of the G / breaking (for example, in A4 it is 
obtained by allowing only some transformation properties for the flavons with non vanishing VEV's). The 
explicit example of A4 is discussed in the next section. Note that, along the same line, a model with fi — r 
symmetry can be realized in terms of the group S3 generated by products of A23 and T (see, for example, 
Ref. 



3 The A 4 Group 

A4 is the group of the even permutations of 4 objects. It has 4!/2=12 elements. Geometrically, it can 
be seen as the invariance group of a tetrahedron (the odd permutations, for example the exchange of two 
vertices, cannot be obtained by moving a rigid solid). Let us denote a generic permutation (1, 2, 3, 4) — > 
(ni, 7i2, ?i3, simply by (12^211-^4) . A4 can be generated by two basic permutations S and T given by 
S = (4321) and T = (2314). One checks immediately that 

S 2 = T 3 = (ST) 3 = 1 . (16) 

This is called a "presentation" of the group. The 12 even permutations belong to 4 equivalence classes (h 
and k belong to the same class if there is a g in the group such that ghg^ 1 = k) and are generated from S 
and T as follows: 

Ci : 1= (1234) 

C 2 : T = (2314), ST = (4132), TS = (3241), STS = (1423) 

(11) 

C 3 : T 2 = (3124), ST 2 = (4213), T 2 S = (2431), TST = (1342) 
C 4 : S = (4321), T 2 ST = (3412), TST 2 = (2143) 

Note that, except for the identity I which always forms an equivalence class in itself, the other classes are 
according to the powers of T (in C4, S could as well be seen as ST 3 ). 
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Table 3 Characters of A4 



Class 


X 1 


x 1 ' 


1" 

X 


X 3 


Ci 


1 


1 


1 


3 


c 2 


1 




UJ 2 





c 3 


1 


w 2 









1 


1 


1 


-1 



The characters of a group are defined, for each element g, as the trace of the matrix that maps 
the element in a given representation R. From the invariance of traces under similarity transformations it 
follows that equivalent representations have the same characters and that characters have the same value 
for all elements in an equivalence class. Characters satisfy ^2 g Xg^Xg* — N5 RS ', where N is the number 
of transformations in the group (N = 12 in A4). Also, for each element h, the character of h in a direct 
product of representations is the product of the characters: = XhX'h an d ^ so i s e q ua l to the sum 

of the characters in each representation that appears in the decomposition of R (8 S. In a finite group the 
squared dimensions of the inequivalent irreducible representations add up to N. The character table of 
A4 is given in Tab. [3] From this table one derives that A 4 has four inequivalent representations: three of 
dimension one, 1,1' and 1" and one of dimension 3. 

It is immediate to see that the one-dimensional unitary representations are obtained by: 

1 S =1 T = 1 

I' S=l T = e i27T / S = oj (18) 

1" S=l T = e iAlr / 3 = cj 2 . 

Note that u = — 1/2 + iy/3/2 is the cubic root of 1 and satisfies uj 2 = lo*, 1 + uj + uj 2 = 0. 

The three-dimensional unitary representation, in a basis where the element S — S' is diagonal, is built 
up from: 



S' = -1 , T' = 1 . (19) 





The multiplication rules are as follows: the product of two 3 gives 3x3 = 1 + V + 1" + 3 + 3 
and V X V = 1", V x 1" = 1, 1" x 1" = V etc. If 3 — (ai,a 2 ,a 3 ) is a triplet transforming by 



the matrices in eq. ( 19i we have that under S'; S'(ai, 0,2, 03)' = (ai, — 0,2, —03)' (here the upper index 
t indicates transposition) and under T': T'(ai, aa, 013)* = (02,03,01)*. Then, from two such triplets 
3 a ~ (a\, a 2 , 03), 3f, ~ (pi, 62, &s) the irreducible representations obtained from their product are: 

1 = ai&i + a 2 b 2 + 0363 

l' = CJ1&1 + W 2 a2&2 + wa 3^3 

1" = a\bi + wa 2 &2 + w 2 a 3 6 3 (20) 
3 ~ (a2b 3 ,a 3 b 1 ,a 1 b 2 ) 
3 ~ (a 3 b2,aib 3 ,a 2 bi) 

In fact, take for example the expression for 1" = a\b\ +wa2&2 +w 2 a3&3. Under S' it is invariant and under 
T' it goes into 0262 + ^03^3 + co 2 aibi = u 2 [a\bi + ^02^2 + u 2 a 3 b 3 ] which is exactly the transformation 
corresponding to 1". 
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In eq. ( 19 1 we have the representation 3 in a basis where S is diagonal. We shall see that for our purposes 
it is convenient to go to a basis where instead it is T that is diagonal. This is obtained through the unitary 
transformation: 



T = VT'V^ = 



S 



VS'V^ = 




(21) 



(22) 



where: 



V 




(23) 



The matrix V is special in that it is a 3 x 3 unitary matrix with all entries of unit absolute value. It is 
interesting that this matrix was proposed long ago as a possible mixing matrix for neutrinos 187,88]. We 
shall see in the following that in the T diagonal basis the charged lepton mass matrix (to be precise the 
matrix m\m e ) is diagonal. Notice that the matrices (S, T) of eqs. (21 1 and (22 1 coincide with the matrices 
(Stb, Ttb) of the previous section. 

In this basis the product rules of two triplets, (ipi,^,^) an d 1,^2,^3) of A> according to the 
multiplication rule 3x3 = 1 + 1' + 1" + 3 + 3 are different than in the S diagonal basis (because for 
Majorana mass matrices the relevant scalar product is (ab) and not (a^6))and are given by: 



- 



■ + i>2<Pl ' 



V>2^2 + ?/> 3 ^l + i>lVZ 



2^1 </Ji - t/>2</?3 

2i/W3 - ^1^2 




1, 




1', 




1", 




^2^3 


- ^3^2 


VW2 


- 1p2<fl 




- Vl^3 



(24) 



3a 



An obvious representation of A 4 is obtained by considering the 4x4 matrices that directly realize each 
permutation. For S = (4321) and T = (2314) we have 



S 1 



/ 1 \ 

10 

10 

y 1 j 



T 4 



/ 





1 








\ 








1 









1 













V 











1 





(25) 



The matrices S4 and T4 satisfy the relations in eq. ( 16 1, thus providing a representation of At. Since the 



only irreducible representations of A are a triplet and three singlets, the 4x4 representation described by 
£4 and T4 is not irreducible. It decomposes into the sum of the invariant singlet plus the triplet representa- 
tion. In fact the vector (1,1,1, 1) T is clearly invariant under permutations and similarly the 3-dimensional 
space orthogonal to it. In matrix terms this decomposition is realized by the unitary matrix [89] U given 
by 



U =2 



f +1 

-1 

+ 1 

V +1 



+1 +1 \ 

+1 -1 

+1 -1 

-1 -1 / 



(26) 
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This matrix maps S4 and T4 into matrices that are block-diagonal: 



/ JJ \ 



us 4 u^ = 



s 



V 



(A 



UT 4 U^ = 



) 



\ 



T 

J 



(27) 



where S and T are the generators of the three-dimensional representation in eq. ( 19 1. 

In the following we will work in the T diagonal basis, unless otherwise stated. In this basis the 12 
matrices of the 3-dimensional representation of A4 are given as follows: 



d: 1 = 



T 



TS = 



-3 • 




2 

2uj 



--4 • 



-1 


2 


2 


2a; 2 


-u? 


2uj 2 


2uj 


2tu 


— UJ 


1 2 


2 




-1 


2 


)' 


2 


-1 




/ 1 


2oj 2 


2lo 


2w 


-1 


2uj 2 


V 2uj 2 


2uj 


-1 



ST 



STS = 



ST 2 



TST = 



T 2 ST = - 






2a; 2 


2a; 


2a; 2 


—a; 




2a; 


2 


-a; 2 


-1 


2a; 2 


2a; \ 


2 


-UJ 2 


2a; 


2 


2a; 2 


-u J 




2a; 


2uj 2 


2cj 


-w 2 


2 


x 2w 2 


2 


—UJ 




2w 


2uj 2 


f 2J 2 


-1 


2uj 


\ 2a; 


2a; 2 


-1 



T 2 



T 2 S = 



s- 1 



TST 2 = 



We can now see why A4 works for TB mixing. In Sec. [2] we have already mentioned that the most gen- 
eral mass matrix for TB mixing in eq. ([8]), in the basis where charged leptons are diagonal, can be specified 
as one which is invariant under the 2-3 (or /i — t) symmetry and under the S unitary transformation, as 
stated in eq. ( 12 1. This observation plays a key role in leading to A4 as a candidate group for TB mixing, 
because S is a matrix of A4. Instead the matrix A23 is not an element of A4 (because the 2-3 exchange is 
an odd permutation). We shall see that in A4 models the 2-3 symmetry of the neutrino mass matrix arises 
as an accidental symmetry of the LO Lagrangian by imposing that there are no flavons transforming as 1' 
or 1" that break A 4 with two different VEV's (in particular one can assume that there are no flavons in the 
model transforming as 1' or 1"). It is also clear that a generic diagonal charged lepton matrix m\m e is 
characterized by the invariance under T, or T^m\m e T — m\m e . 

The group A4 has two obvious subgroups: Gs, which is a reflection subgroup generated by S, and Gt, 
which is the group generated by T, which is isomorphic to Z3. If the flavour symmetry associated to A4 is 
broken by the VEV of a triplet ip = (ipi , ip2, ^3) of scalar fields, there are two interesting breaking pattern. 
The VEV 



(ip) = (v s ,v s ,v s ) 



(28) 
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breaks A4 down to Gs, while 

(<p) = (v T ,0,Q) (29) 

breaks A4 down to Gt- As we will see, Gs and Gt are the relevant low-energy symmetries of the neutrino 
and the charged-lepton sectors, respectively. Indeed we have already seen that the TB mass matrix is 
invariant under Gs and a diagonal charged lepton mass m\m e is invariant under Gt- 

4 Applying A 4 to Lepton Masses and Mixings 



In the lepton sector a typical A4 model works as follows [90) . One assigns leptons to the four inequivalent 
representations of A 4 : LH lepton doublets I transform as a triplet 3, while the RH charged leptons e c , 
H c and r c transform as 1, 1" and V, respectively. These models can be realized both with and without a 
see-saw mechanism. In the first case there are three right-handed neutrinos transforming as a triplet of A4, 
while in the second case the source of neutrino masses is a set of higher dimensional operators violating 
the total lepton number. Here we consider a see-saw realization, so we also introduce conjugate neutrino 
fields v c transforming as a triplet of A 4 . The fact that LH lepton doublets I and, in the see-saw case, also 
the RH neutrinos v c , transform as triplets is crucial to realize the fixed ratios of mass matrix elements 
needed to obtain TB mixing. A drawback is that for the ratio r, defined by A?rig Un /Am^ tm , one would 
expect y/r w O(l) to be compared with the experimental value is y/r « 0.2, which implies a moderate 
fine-tuning. 

One adopts a supersymmetric (SUSY) context also to make contact with Grand Unification (flavour 
symmetries are supposed to act near the GUT scaleQ. In fact, as well known, SUSY is important in 
GUT's for offering a solution to the hierarchy problem, for improving coupling unification and for making 
the theory compatible with bounds on proton decay. But, in models of lepton mixing, SUSY also helps 
for obtaining the vacuum alignment, because the SUSY constraints are very strong and limit the form of 
the superpotential very much. Thus SUSY is not necessary but it is a plausible and useful ingredient. The 
flavour symmetry is broken by two sets of flavons $ e and <!>„, invariant under the SM gauge symmetry, that 
at the LO break A4 down to Gt and Gs, respectively. At this order $ e couples only to the charge lepton 
sector and to the neutrino sector. Typically <l> e and <!>„ include triplets and invariant singlets under A4, 
but models with flavons transforming as 1' and 1" have also been considered [96,97]. For example $ e can 
consist of the triplet (pT with the vacuum alignment in eq. (|29]i and can include the triplet ips with the 



vacuum alignment in eq. (28 1 and two invariant singlets £, £. Two Higgs doublets h Ut d, invariant under 
A4, are also introduced. One can obtain the observed hierarchy among m e , m M and m r by introducing 
an additional U(l)^ n flavour symmetry J98) under which only the RH lepton sector is charged (recently 
some models were proposed with a different VEV alignment such that the charged lepton hierarchies are 
obtained without introducing a U(l) symmetry | |99||100) ). We recall that U(1)fn is a simple flavour 
symmetry where particles in different generations are assigned (in general) different values of an Abelian 
charge. Also Higgs fields may get a non zero charge. When the symmetry is spontaneously broken the 
entries of mass matrices are suppressed if there is a charge mismatch and more so if the corresponding 
mismatch is larger. We assign FN-charges 0, q and 2q to t c , /i c and e c , respectively. There is some 
freedom in the choice of q. Here we take q — 2. By assuming that a flavon 9, carrying a negative unit of 
FN charge, acquires a VEV (6>)/A = Ac < 1, where Ac = sin 6c, the Yukawa couplings become field 
dependent quantities y e . M . T = y e . M . r (^) and we have 

lfr«0(l), y M «0(A c ), y e «0(A c ). (30) 

Had we chosen q = 1, we would have needed (0)/A of order X c , to reproduce the above result. The 
superpotential term for lepton masses, wi is given by 

wi = y e e c (^ T 0+y M Ai c (^T0'+yrT c (^T0''+y(^ c + (^e+5A|)(^ c ^ c )+^i3(^s^ c ) + - ■ ■ (31) 



^ When the flavour symmetry is broken contextually with the electroweak one, such as in Refs. |9l]}93|, strong constraints from 



FCNC transitions are usually present 1 94 95], that can eventually rule out the model. 
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with dots denoting higher dimensional operators that lead to corrections to the LO approximation. In our 
notation, the product of 2 triplets (33) transforms as 1, (33)' transforms as 1' and (33)" transforms as 
1". To keep our formulae compact, we omit to write the Higgs and flavon fields h Uj d, & and the cut-off 
scale A. For instance y e e c (ipTl) stands for y e e c (ipTl)hd0 4: / A 5 . The parameters of the superpotential Wi 
are complex, in particular those responsible for the heavy neutrino Majorana masses, xj^ q. Some terms 
allowed by the A4 symmetry, such as the terms obtained by the exchange (fx *H> ips, (or the term (u c u c )) 
are missing in wi. Their absence is crucial and, in each version of A4 models, is motivated by additional 
symmetries. 

The LO superpotential in eq. |3lj ) leads to a diagonal mass matrix for the charged lepton^] 




m e 0) =v d I Vlx J n with rj = ~ , (32) 

and to a neutrino mass matrix mi?' of the same form as that of eq. (Ji}. As for the neutrino spectrum both 
normal and inverted hierarchies can be realized. It is interesting that A4 models with the see-saw mech- 
anism typically lead to a light neutrino spectrum which satisfies the sum rule (among complex masses): 



— = - '-■ (33) 

m.2, 7Yl\ TO 2 



A detailed discussion of a spectrum of this type can be found in Refs. 1 90 100- 102|. The above sum rule 



gives rise to bounds on the lightest neutrino mass. As a consequence, for example, the possible values of 
\m ee \ are restricted. For normal hierarchy we have 

Am 2 sun w 0.007 eV . (34) 



1 ee| 3^3 
while for inverted hierarchy 



Am ltm 



0.017 eV. (35) 



In a completely general framework, without the restrictions imposed by the flavour symmetry, |m ee | could 
vanish in the case of normal hierarchy. In this model \m ee \ is always different from zero, though its value 
for normal hierarchy is probably too small to be detected in the next generation of 0v(3(3 experiments. 

In the leading approximation A4 models lead to exact TB mixing. In these models TB mixing is implied 
by the symmetry at the leading order approximation which is corrected by non-leading effects. Given the 
set of flavour symmetries and having specified the field content, the non leading corrections to TB mixing, 
arising from higher dimensional effective operators, can be evaluated in a well defined expansion. 

The departure from the LO approximation depends on the subleading contributions Srrdp, £mi , to 
the charged lepton and the neutrino mass matrices, respectively: 

m e = m ( 0) + 8m& + ... , m v = + Smfi + ... , (36) 



which can vary according to the model considered. In all models considered here [59 90 100 103 1 the 
NLO corrections to the charged lepton mass matrix are of the following type: 

0(y e ) 0(y e 
SmP=v d ( 0(y M ) 0(y„) O(y^) | ^ 
0(y T ) 0(y T ) 




^ We absorbed in yf (f = e,^,,r) the appropriate factor of (6) / A. 
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where £ is small adimensional parameter given by the ratio between the fiavon VEVs and A. The transfor- 
mation needed to diagonalize m e is V^Fm e U e = mf ,a9 where 

U e = -cfU 1 c| 3 £ . (37) 

V -<%z i / 

To discuss the NLO contribution to m v we distinguish two cases. 



4.1 Typical A 4 Models 

In "typical" A 4 models J90|l00|l03 |, the NLO contribution Sm^P in eq. (361) is a generic symmetric matrix 



with entries suppressed, compared to the corresponding entries in rn^ \ by a relative factor of the order 
of the ratio between a fiavon VEV and A. This occurs both with and without the see-saw mechanism. The 
generic transformation that diagonalizes m v is U^U^ B m v UTBUu where 



U v 



1 c^C 

-<fj' f i f 

~ c 13 £ — c 23 £ 1 



(38) 



and C23 are complex parameters of order one in absolute value. Barring a fine-tuning of the 



where c\ 2 

Lagrangian parameters, in these models the suppression factors £ and £' are expected to be of the same 
order of magnitude. For example, beyond the LO the equations satisfied by ($ e ) and ($„) are no longer 
decoupled and the corrections to the LO fiavon VEVs turn out to be of the same size, for both <£> e and $„. 
All the elements of the mixing matrix get corrections of the same size £ « We expecj^On the contrary, 
the expressions in eq. ( |40| show explicitly the dependence of the NLO mixing angles on the corrections 
from both the neutrino and charged lepton sectors.: 



sin 2 6*23 



sin 2 012 = * - \ne{c\ 2 + c e 13 ) £ + ?p ft e (c? 2 ) £ 



(40) 



sm &13 



3\/2(cJ 2 



-13; 



2\/3 f\/2c^ 3 + 



According to these expressions, in order to reach the central value for the reactor angle in agreement 
with eq. the parameter £ is expected to be O(0. 1). A precise value can be found by studying the 
success rate to reproduce all the three mixing angles inside the corresponding 3<r ranges, depending on 



the value of £. As shown in ref. |58 1, in a scan with the cA parameters that multiply £ treated as random 
complex numbers with absolute values following a Gaussian distribution around 1 with variance 0.5, the 
value of £ that maximizes the success rate is found to be 0.075(0.078) for the NH (IH). The corresponding 
success rate is ~ 12%, which is not large but not hopelessly small either. For this value of £ in Fig. [2] we 
quantitatively analyze the expressions in eq. (pOl and their correlations: in the plots on the left (right), we 
show the correlation between sin 2 #i 3 and siir6i 2 (sin 2 6*23). In the plots we show only the NH case. The 
IH case is similar. 

As we can see, the plots are representing the general behaviour of this class of models: sin 2 (9 13 increases 
with £, but correspondingly also the deviation of sin 2 9i 2 from 1/3 does. As a result, even for the value 



Eq. |40| is a particular case of the general parametrization presented Ref. [104]: 

1 1 r 

sin0 2 3 = — t=(1 + a ) ) sin0i2 = — ^=(1 + s) , sin 6*13 ; 



(39) 



with a, s and r real numbers. The expressions in Eq. (40) show explicitly the dependence of the NLO mixing angles on the corrections 
from both the neutrino and the charged lepton sectors. 
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Fig. 2 Typical A4 Models. On the left (right), we plot sin 2 #13 as a function of sin 2 #12 (sin 2 #23), following 
eq. \4Q\ . The dashed-black lines represent the 3<r values for the mixing angles from the Fogli et al. fit [7|. Only the 
NH data sets is shown. The parameter £ is taken equal to 0.075. The Cy" parameters that multiply £ are treated as 
random complex numbers with absolute values following a Gaussian distribution around 1 with variance 0.5. 



of £ that maximizes the success rate, the requirement for having a reactor angle inside its 3<r error range 
corresponds to a prediction for the solar angle that spans all the 3cr experimental error bar and is often not 
even in agreement with the data. 



4.2 Special A 4 Models 



In these models Srriv in eq. ( 36 1 is protected by the assumed symmetries so that it remains invariant under 



S and its relative size, compared to can be bigger than £. For instance, in Ref. [ |59[ (<J> e ) an d {&u) 
satisfy decoupled equations up to NLO so that it is possible to achieve ($ e ) < ($„). Moreover ($„) 
couples to charged leptons only at the NNLO so that the dominant source of corrections to the neutrino 
mixing pattern is 8m\}\ Including these NLO corrections we have 

x y — w y + w \ 

y — w x + z + w y — z , (41) 
y+w y—z x+z—w J 

where w represents the part of the NLO corrections that cannot be absorbed by a redefinition of x, y and z. 
The parameter w is smaller that x, y, z, but not necessarily much smaller. The crucial property of 77i„ |58| 
is that it is still invariant under the action of S (but not any more of A23): 

S T m v S — m v . (42) 

Actually it can be directly proven that the matrix in eq. ( |4T| ) is the most general one invariant under S. The 
matrix m v can be diagonalized in two steps. First we transform m v by a Tri-Bimaximal rotation: 

(x — y \/3w \ 

x + 2y , (43) 

y/3w x-y + 2z J 

Second, we perform a unitary transformation in the (1,3) plane: 

m 2 + ici 2 = i, (44) 
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V T m' u V = m d ; ag 
The exact rotation is given by: 

2a£' uv*(u* 



|a| 2 -|e| 2 \v\ 2 -\u\ 2 ' 
The unitary matrix that diagonalizes m v is 

U TB V = 



with u 



2^/3w 
x-y 



and v 



x — y + 2z 



(45) 



(46) 



^J2j3a 1/V3 
-a/V6 + £*/V2 1/V3 -a*/V2-C/V6 
-a/V6-£*/V2 1/V3 +a*/y/2-S'/y/6 



(47) 



It is not restrictive to choose a real and positive and we have: In eq. ( |47| it is not restrictive to choose a 
real and positive and we have: 



S C p ~ 

sin 6 13 = 

sin 2 6>i2 = 

sin 2 6*23 = 



arg£' 
2 



3-2|£' 



V2 
3' 

1 



jjia 2 -^e(c? 2 + C f 3 K 



(48) 
(49) 

(50) 



V3a 



1 



73 



ir|cos5 + 7ee(c| 3 )e (51) 



where we have also included the effects coming from the diagonalization of the charged lepton sector, to 
first order in £. The second equality shows the result expanded in powers of |£'|, to the order |£'| 2 . In these 



models is of order 0.1, bigger than £ so that the contribution of eq. (37 i are subdominant 



It is interesting to note that if we neglect the corrections proportional to £, we have an exact relation 
between the solar and the reactor anglaj 



sin 2 0i2 



1 



3(1 - sin^ 0i 3 



1 



V2 



sin 0i3 cos Sep 



(52) 



The first expression demonstrates that the unitary transformation V always increases the solar angle from 
the TB value, while the preferred la interval is below the TB prediction. This is a small effect, of second 
order in 0i 3 , that can be compensated by the corrections proportional to £. The second correlation involves 
the Dirac CP phase and is particularly interesting considering the recent hint of a CP phase close to tt 
for the NH case: when considering the la (2a) ranges for the mixing angles, one sees an indication that 
cos<5cp lies m the interval [—1, —0.5], while no indication arises when the 3a error band for sin 2 023 is 
taken into account. Although these results for the CP phase is modified by the inclusion of the subleading 
£ contributions, these correlations will allow an interesting test for such models once Sep is measured and 
the precision on sin 2 023 is improved. 

The success rate to reproduce all the three mixing angles inside their corresponding 3a error ranges, 
as a function of |£'|, is studied in ref. ]58) . The parameters are chosen such that £ is a real number in 
[0.005, 0.06] and are random complex numbers with absolute values following a Gaussian distribution 
around 1 with variance 0.5. The value of that maximizes the success rate for both the hierarchies is 



It has been shown in Ref. [ 105], from general group theoretical considerations, that these correlations are a general feature of 
flavour models when the symmetry group of the charged lepton (neutrino) mass matrix is Z3 (Z2). 
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found to be 0.183. The corresponding success rate is much larger in these models (~ 64%) than for the 



typical A4 models. For the stated range of £ and the optimal value of £' the deviations in eqs. (50 1 and (51 



and their correlations are quantitatively analyzed in Fig. [3] in the plots on the left (right) column, we show 



the correlations in eqs. (50i and (51 1 between sin #13 and sin 812 or sin 823, respectively. We see that, 
for this choice of the parameters, the model can well describe all three angles inside the corresponding 3<7 
interval, and its success rate is much larger than that of the typical TB models. 




Fig. 3 Special A4 Models, sin 2 #13 as a function of sin 2 612 (sin 2 623 is plotted on the left (right), following 
eqs. l |50[ > and (T5TJ. The dashed-black lines represent the 3<r values for the mixing angles from the Fogli et al. fit (7J. 
Only the NH data sets is shown. The parameter £ is a real number in [0.005, 0.06]; £' is a complex number with 
absolute values equal to 0.183; the parameters Cy are random complex numbers with absolute values following a 
Gaussian distribution around 1 with variance 0.5. 



One particularly interesting example realizing this scenario is provided by the Lin model [59] (see 
also JTOrj) ), formulated before the T2K, MINOS, DOUBLE CHOOZ, Daya Bay and RENO results were 
known. In the Lin model, the A4 symmetry breaking is arranged, by suitable additional Z n parities, in 
such a way that, not only at LO but also at NLO, the corrections to the charged lepton and the neutrino 
sectors are kept separate. Then the contributions to neutrino mixing from the diagonalization of the charged 
leptons can be of O(Xq) while those in the neutrino sector can be of 0(Xc)- In addition, in the Lin model 
these large corrections do not affect 6*12 and satisfy the relations in eq. pi) , with 6 being the CKM-like 
CP violating phase of the lepton sector. Thus in the Lin model the NLO corrections to the solar angle 6*12 
and to the reactor angle 6*13 can naturally be of different orders and 6*13 ~ O(Xc) is not at all surprising. 

A related scenario is provided by a framework based on a typical A4 model as described at the beginning 
of this section, extended by the inclusion in <!>„ of two additional singlets transforming as V and 1" [ 107 1. 
Once these singlets and the triplet tpg acquire a VEV, the resulting LO neutrino mass matrix is only 
invariant under the action of S and no more under A23. Thus, already at the LO we have of the 
general form in eq. (41 1. In this framework the smallness of #13 is however unexplained. 



4.3 Summary on A4 models for lepton mixing 

In summary, in the absence of specific dynamical tricks, in a generic A4 model, all the three mixing angles 
receive corrections of the same order of magnitude. Since the experimentally allowed departures of 612 
from the TB value, sin 2 612 = 1/3, are small, numerically not larger than O(Xq), it follows that both 
013 and the deviation of #23 from the maximal value are expected to also be typically of the same general 
size. The central values sin #13 f=s 0.15 that can be derived from the experimental results in Tab. [2] are in 
between 0{X 2 C ) ~ O(0.05) and 0(X C ) ~ 0(0.23). Although models based on TB (or GR) mixing tend 
to lead to a smaller value of #13 one can argue that they are still viable with preference for the lower side 
of the experimental range. But, as we have seen, one can introduce some additional theoretical input to 
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enhance the value of 6*13 in an A4 model (several models have been recently proposed in order to fulfill 
this goal [86 97||105]|107 ■ 153| ). As a result we now have examples of A4 models where the departures 
from exact TB mixing are naturally larger for 6*13 than for #12 and 6*23. 



5 A 4 , Quarks and GUT's 

Much attention has been devoted to the question whether models with TB mixing in the neutrino sector can 
be suitably extended to also successfully describe the observed pattern of quark mixings and masses and 
whether this more complete framework can be made compatible with (supersymmetric) SU(5) or SO(10) 
Grand Unification. 

The simplest attempts of directly extending models based on A4 to quarks have not been satisfactory. 
At first sight the most appealing possibility is to adopt for quarks the same classification scheme under A4 
that one has used for leptons (see, for example, Ref. [90 1). Thus one tentatively assumes that LH quark 
doublets Q transform as a triplet 3, while the antiquarks (u c ,d c ), (c c , s c ) and (t c ,b c ) transform as 1, 1" 
and 1', respectively. This leads to V u = Vd and to the identity matrix for Vckm — V^Vd in the lowest 
approximation. This at first appears as very promising: a LO approximation where neutrino mixing is 
TB and Vckm = 1 is a very good starting point. But there are some problems. First, the corrections to 
Vckm = 1 turn out to be strongly constrained by the leptonic sector, because lepton mixing angles are 
very close to the TB values, and, in the simplest models, this constraint leads to a too small V us (i.e. the 
Cabibbo angle is rather large in comparison to the allowed shifts from the TB mixing angles). Also in 
these models, the quark classification which leads to Vckm = 1 is not compatible with A4 commuting 
with SU(5). An additional consequence of the above assignment is that the top quark mass arises from a 
non-renormalizable dimension-5 operator. In that case, to reproduce the top mass, we need to compensate 
the cutoff suppression by some extra dynamical mechanism. Alternatively, we have to introduce a separate 
symmetry breaking parameter for the quark sector, sufficiently close to the cutoff scale. 

Due to this, larger discrete groups have been considered for the description of quarks. A particularly 
appealing set of models is based on the discrete group T', the double covering group of A4 | 154 162 1. 



The representations of T' are those of A4 plus three independent doublets 2, 2' and 2". The doublets 
are interesting for the classification of the first two generations of quarks [ 163f [165) . For example, in 
Ref. |155| a viable description was obtained, i.e. in the leptonic sector the predictions of the A4 model 
are maintained, while the T' symmetry plays an essential role for reproducing the pattern of quark mixing. 
But, again, the classification adopted in this model is not compatible with Grand Unification. 

As a result, the group A4 was considered by many authors to be too limited to also describe quarks and to 
lead to a grand unified description. But it has been shown [ 166 1 that this negative attitude is not justified and 
that it is actually possible to construct a viable model based on A4 which leads to a grand unified theory 
(GUT) of quarks and leptons with TB mixing for leptons and with quark (and charged lepton) masses 
and mixings compatible with experiment. At the same time this model offers an example of an extra 
dimensional SU(5) GUT in which a description of all fermion masses and mixings is accomplished. The 
formulation of SU(5) in extra dimensions has the usual advantages of avoiding large Higgs representations 
to break SU(5) and of solving the doublet-triplet splitting problem. The choice of the transformation 
properties of the two Higgses H5 and has a special role in this model. They are chosen to transform 
as two different A 4 singlets 1 and 1'. As a consequence, mass terms for the Higgs colour triplets are 
not directly allowed and their masses are introduced by orbifolding, a la Kawamura |167| . In this model, 
proton decay is dominated by gauge vector boson exchange giving rise to dimension-6 operators, while 
the usual contribution of dimension-5 operators is forbidden by the selection rules of the model. Given 
the large Mcut scale of SUSY models and the relatively huge theoretical uncertainties, the decay rate is 
within the present experimental limits. A see-saw realization in terms of an A 4 triplet of RH neutrinos 
v c ensures the correct ratio of light neutrino masses with respect to the GUT scale. In this model extra 
dimensional effects directly contribute to determine the flavour pattern, in that the two lightest tenplets T\ 
and T2 are in the bulk (with a doubling Ti and T[, i = 1, 2 to ensure the correct zero mode spectrum), 
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whereas the pentaplets F and T3 are on the brane. The hierarchy of quark and charged lepton masses and 
of quark mixings is determined by a combination of extra dimensional suppression factors and of U{\)fn 
charges, both of which only apply to the first two generations, while the neutrino mixing angles derive from 
A4 in the usual way. If the extra dimensional suppression factors and the U(1)fn charges are switched 
off, only the third generation masses of quarks and charged leptons survive. Thus the charged fermion 
mass matrices are nearly empty in this limit (not much of A4 effects remain) and the quark mixing angles 
are determined by the small corrections induced by the above effects. The model is natural, since most 
of the small parameters in the observed pattern of masses and mixings as well as the necessary vacuum 
alignment are justified by the symmetries of the model. However, in this case, like in all models based 
on U(1)fn> the number of 0(1) parameters is larger than the number of measurable quantities, so that 
in the quark sector the model can only account for the orders of magnitude (measured in terms of powers 
of an expansion parameter) and not for the exact values of mass ratios and mixing angles. A moderate 
fine-tuning is only needed to enhance the Cabibbo mixing angle between the first two generations, which 
would generically be of O(Xq). 

The problem of constructing GUT models based on SU(5) <g> Gf or 5*0(10) ® Gt with approximate TB 
mixing in the leptonic sector has also been considered by many authors. Examples are: for Gf = A4 Ref. 
| 166l[T68||T75) , for V Ref. fl5^[l6l) , for 5 4 Ref. ||176ffl79|. A s for the models based on 50(10) ® G f 



recent examples were discussed with G/ = 5 4 | |54||146||180||181| and Gf = PSL 2 {7) | |182||183| . Clearly 
the case of 50(10) is even more difficult than that of SU (5) because the neutrino sector is tightly related 
to that of quarks and charged leptons as all belong to the 16 of 50(10). For a discussion of 50(10) ® A4 
models, see |184| . More in general see Refs. [119 185— 188]. In our opinion most of the models are 
incomplete (for example, the crucial issue of VEV alignment is not really treated in depth as it should) 
and/or involve a number of unjustified steps and ad-hoc fine-tuning of parameters. 



6 Possible Origin of A 4 

There is an interesting relation [90 1 between the A4 model considered so far and the modular group. This 
relation could possibly be relevant to understand the origin of the A4 symmetry from a more fundamental 
layer of the theory. The modular group T is the group of linear fractional transformations acting on a 
complex variable z: 

z — > — , ad — be — 1 , (53) 

cz + a 

where a, b, c, d are integers. There are infinite elements in T, but all of them can be generated by the two 
transformations: 

s : z -> — , t : z -> z + 1 , (54) 

z 



The transformations s and t in ( 54 1 satisfy the relations 

s 2 = {stf = 1 (55) 
and, conversely, these relations provide an abstract characterization of the modular group. Since the re- 



lations in eqs. ( 16 1 are a particular case of the more general constraint in eq. (55 1, it is clear that A4 is 
a very small subgroup of the modular group and that the A4 representations discussed above are also 
representations of the modular group. In string theory the transformations in eq. |54]i operate in many 
different contexts. For instance the role of the complex variable z can be played by a field, whose VEV 
can be related to a physical quantity like a compactification radius or a coupling constant. In that case s in 



eq. ( 54 1 represents a duality transformation and t in eq. p4\ represents the transformation associated to an 



axionic symmetry. 
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A different way to understand the dynamical origin of A4 was presented in Ref. [89| where it is shown 
that the A4 symmetry can be simply obtained by orbifolding starting from a model in 6 dimensions (6D). 
In this approach A4 appears as the remnant of the reduction from 6D to 4D space-time symmetry induced 
by the special orbifolding adopted. There are 4D branes at the four fixed points of the orbifolding and 
the tetrahedral symmetry of A4 connects these branes. The standard model fields have components on the 
fixed point branes while the scalar fields necessary for the A4 breaking are in the bulk. Each brane field, 
either a triplet or a singlet, has components on all of the four fixed points (in particular all components are 
equal for a singlet) but the interactions are local, i.e. all vertices involve products of field components at 
the same space-time point. This approach suggests a deep relation between flavour symmetry in 4D and 
space-time symmetry in extra dimensions. 

The orbifolding is defined as follows. We consider a quantum field theory in 6 dimensions, with two 
extra dimensions compactified on an orbifold T 2 JZ^. We denote by z = x§ + i x§ the complex coordinate 
describing the extra space. The torus T 2 is defined by identifying in the complex plane the points related 
by 

z ->• z + 1 

(56) 



■7 l = e^'\ 



where our length unit, 2ttR, has been set to 1 for the time being. The parity Z2 is defined by 

z -> -z (57) 

and the orbifold T 2 jZ% can be represented by the fundamental region given by the triangle with vertices 
0,1,7, see Fig. |4] The orbifold has four fixed points, {z\, z 2 , z 3 , z 4 ) = (1/2, (1 + 7)/2, 7/2, 0). The 
fixed point Z4 is also represented by the vertices 1 and 7. In the orbifold, the segments labelled by a 
in Fig. |4] (0, 1/2) and (1, 1/2), are identified and similarly for those labelled by 6, (1, (1 + j)/2) and 
(7, (1 + 7) /2), and those labelled by c, (0, 7/2), (7, 7/2). Therefore the orbifold is a regular tetrahedron 
with vertices at the four fixed points. 




Fig. 4 Orbifold T2 /Z2 ■ The regions with the same numbers are identified with each other. The four triangles bounded 
by solid lines form the fundamental region, where also the edges with the same letters are identified. The orbifold 
T2 /Z2 is exactly a regular tetrahedron with 6 edges a, b, c, d, e, f and four vertices Zi, 22, £3, 24, corresponding to the 
four fixed points of the orbifold. 



The symmetry of the uncompactified 6D space time is broken by compactification. Here we assume 
that, before compactification, the space-time symmetry coincides with the product of 6D translations and 
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6D proper Lorentz transformations. The compactification breaks part of this symmetry. However, due to 
the special geometry of our orbifold, a discrete subgroup of rotations and translations in the extra space is 
left unbroken. This group can be generated by two transformations: 

S :z z + - 

2 (58) 

T :z — >• ujz w = 7 2 . 

Indeed S and T induce even permutations of the four fixed points: 

S :{zi,z 2 , z 3 , z±) -> (z 4 , z 3 , z 2l zi) 
T :(zi, z 2 , z 3 , za) -t [z 2 , z 3 , zi, Zi) , 

thus generating the group Aj. From the previous equations we immediately verify that S and T satisfy 
the characteristic relations obeyed by the generators of A4. S 2 = T 3 = (ST) 3, = 1. These relations are 
actually satisfied not only at the fixed points, but on the whole orbifold, as can be easily checked from the 



general definitions of S and T in eq. (58 1, with the help of the orbifold defining rules in eqs. (56 1 and ( |57| i. 
We can exploit this particular geometry of the internal space to build a model with A4 flavour symmetry. 
There are 4D branes at the four fixed points of the orbifolding and the tetrahedral symmetry of A4 connects 
these branes. The standard model fields have components on the fixed point branes while the scalar fields 
necessary for the A 4 breaking are in the bulk. Each brane field, either a triplet or a singlet, has components 
on all of the four fixed points (in particular all components are equal for a singlet) but the interactions are 
local, i.e. all vertices involve products of field components at the same space-time point. In the low-energy 
limit this model coincides with one of those presented in sec. |4] [ 90 1 . Unfortunately in such a limit the 6D 
construction does not provide additional constraints or predictions. 

This construction can be embedded in a SU(5) GUT [189) . Other discrete groups can arise from 
the compactification of two extra dimensions on orbifolds and the possibilities have been classified in 
Ref. 1 190 1 within a field theory approach. In string theory the flavour symmetry can be larger than the 
isometry of the compact space. For instance in heterotic orbifold models the orbifold geometry combines 
with the space group selection rules of the string, as shown in Ref. ]191| . Discrete flavour symmetries from 
magnetized/intersecting D-branes are discussed in Ref. ||192|. Discrete symmetries can also arise from the 



spontaneous breaking of continuous ones. Such apossibility has been discussed in Refs. [193 195). 



7 Alternative routes to TB mixing 

While A4 is the minimal flavour group leading to TB mixing, alternative flavour groups have been studied 
in the literature and can lead to interesting variants with some specific features. 

In Ref. |196) , the claim was made that, in order to obtain the TB mixing "without fine-tuning", the finite 
group must be S4 or a larger group containing S4. For us this claim is not well grounded being based on 
an abstract mathematical criterium for a natural model (see also Ref. | |197) ). For us a physical field theory 
model is natural if the interesting results are obtained from the most general lagrangian compatible with 
the stated symmetry and the specified representation content for the flavons. For example, in Ref. 1 90 103 1, 
a natural (in our sense) model for the TB mixing is built with A 4 (which is a subgroup of 54) by simply 
not including symmetry breaking flavons transforming like the 1' and the 1" representations of A4. This 
limitation on the transformation properties of the flavons is not allowed by the rules specified in Ref. [ 196 1, 



which demands that the symmetry breaking is induced by all possible kinds of flavons (note that, according 
to this criterium, the SM of electroweak interactions would not be natural because only Higgs doublets are 
introduced!). Rather, for naturalness we also require that additional physical properties like the VEV 
alignment or the hierarchy of charged lepton masses also follow from the assumed symmetry and are not 
obtained by fine-tuning parameters: for this actually A4 can be more effective than S4 because it possesses 
three different singlet representations 1,1' and 1". 



Copyright line will be provided by the publisher 



fdp header will be provided by the publisher 



19 



Models of neutrino mixing based on #4 have in fact been studied [ 52 - 5*4] [T0T]| 1 76f [TT8 180, 18T| |198 
208 1 . The group of the permutations of 4 objects S4 has 24 elements and 5 equivalence classes (the 
character table is given in Tab.[4]that correspond to 5 inequivalent irreducible representations, two singlets, 
one doublet, two triplets: li, I2, 2, 3i and 32). Note that the squares of the dimensions of all these 
representations add up to 24. 



Table 4 Characters of S4 



Class 


X(li) 


X(la) 


X(2) 


X(3i) 


X(3 2 ) 




1 


1 


2 


3 


3 


c 2 


1 


1 


2 


-1 


-1 


c 3 


1 


-1 





1 


-1 


c A 


1 


1 


-1 








c 5 


1 


-1 





-1 


1 



For models of TB mixing, one starts from the S4 presentation A 3 = B 4 = (BA 2 ) 2 = 1 and identifies, 
up to a similarity transformation, B 2 = S and A = T, where S and T are given in eqs. ( 13 1 and ( 15 1. In 
this presentation one obtains a realisation of the 3-dimensional representation of S4 where the matrices S 
and A-23 in eq. (12 1, that leave invariant the TB form of m„ in eq. (|8j, as well as the matrix T in eq. ( |15) , 
of invariance for m\m e , all explicitly appear |204|. In S4 the 1' and 1" of A4 are collected in a doublet. 
When the VEV of the doublet flavon is aligned along the Gs preserving direction the resulting couplings 
are 2-3 symmetric as needed. In A4 the 2-3 symmetry is only achieved if the V and 1" VEV's are identical 
(which is the S4 prediction). As discussed in Ref. |204) , in the leptonic sector the main difference between 
A 4 and S4 is that, while in the typical versions of A 4 the most general neutrino mass matrix depends on 
2 complex parameters (related to the couplings of the singlet and triplet fiavons), in S4 it depends on 3 
complex parameters (because the doublet is present in addition to singlet and triplet fiavons). 

An interesting deformation of the TB mixing pattern arises from the series of groups A(6n 2 ) that 
generalize the permutation group S4, isomorphic to A(24). Indeed the LO mixing pattern induced by £4 
is completely determined once the residual symmetries G e and G v in the charged lepton sector and in the 
neutrino sector are specified and the TB mixing corresponds to the choice G e = Z3, G v = Z 2 x Z 2 . By 
adopting as flavour group A(6n 2 ) (n = 4, 
U 
n 



UtbV, with V given in eq. (44 1, \£'/a 



, the same choice of residual symmetries leads instead to 
tan(?r/3n) and no CKM-like CP violation fT20l[l38) . For 
4, 8, the mixing pattern is close to the experimental data. Concrete models based on A(96) can be 
found in ref. [142|. The analysis in Refs. [120, 138] accounts only the cases when the residual symmetry 
in the neutrino sector is G„ = Z 2 x Z 2 ; in Refs jT05|p5|[209) , a more general study has been presented, 
where the residual symmetry is G v — Z 2 , while the second Z 2 component arises accidentally. 

Other flavour groups have been considered for models of TB mixing. Some of them include S4 as a 
subgroup, like PSL 2 (7) (the smallest group with complex triplet representations) 1 182 183 210) , while 
others, like A (27) (which is a discrete subgroup of SU(3)) | |2TTH2T5l or Zi x Z 3 ]216| , have no direct 
relation to S4 (2171 . 

A different approach to TB mixing has been proposed and developed in different versions by S. King 
and collaborators over the last few years p3| 



neutrino mass matrix given in eqs. ( |T0[ > and 
lepton basis: 



11 



217 -220]. The starting point is the decomposition of the 



corresponding to exact TB mixing in the diagonal charged 



(2 



(60) 



j(0, — 1, 1), are the respective columns of Utb 



where d>f = i(2, -1, -1), ^ = ^(1, 1, 1), ^ = ^ 
and rrii are the neutrino mass eigenvalues. Such decomposition is purely kinematical and does not possess 
any dynamical or symmetry content. In the King models the idea is that the three columns of Utb &i are 
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promoted to flavon fields whose VEVs break the family symmetry, with the particular vacuum alignments 
along the directions Eq. (60l directly arises in the see-saw mechanism, m v = m^M^iriD, written 
in the diagonal RH neutrino mass basis, M — diag(Mi, M2, M3) when the Dirac mass matrix is given 
by mj) = (ui^i, where m are mass parameters describing the size of the VEVs. In this 

way, to each RH neutrino eigenvalue Mj, a particular light neutrino mass m, is associated. In the case 
of a strong neutrino hierarchy this idea can be combined with the framework of "Sequential Dominance", 
where the lightest RH neutrino, with its symmetry properties fixes the heaviest light neutrino and so on. 
For no pronounced hierarchy the correspondence between Af, and m, can still hold and one talks of "Form 
Dominance" | |221| . In these models the underlying family symmetry of the Lagrangian G/ is completely 
broken by the combined action of the $3 VEVs, and the flavour symmetry of the neutrino mass matrix 



emerges entirely as an accidental residual symmetry of the quadratic form of eq. ( 60 1 [ 2 1 7 ] . The symmetry 
G f plays a less direct role and the name "Indirect Models" is used by the authors. 

An alternative context in which the TB pattern has been implemented is the Holographic Composite 
Higgs Models: theories in extra dimensions that provide a weakly coupled description of certain 4 di- 
mension composite Higgs models. Neutrinos, as the other fermions and the gauge bosons, represent the 
elementary sector, while the SM Higgs arise from the composite Higgs sector. The connection between 
the two sectors is provided through symmetry-defined couplings. In Refs. [ 222 223], it has been shown 
that the TB mixing can arise in this context thanks to a non-Abelian discrete symmetry, S4, A$, A(96) and 
A(384), acting in both the elementary and composite sectors, broken into certain non-trivial subgroups of 
the original symmetry. 



8 Constraints from lepton flavour violating processes 

As we have discussed in the previous sections the relatively large value of 6*13 introduces a marked de- 
parture from the TB limit, while the values of 812 and 623 are very close to it. One challenge for flavour 
models in the lepton sector is to produce in a natural way a relatively large correction to # 13 without af- 
fecting too much the other mixing angles. Another challenge arises from the existing stringent bounds 
on lepton flavour violating processes. In particular, we refer to the recent improved MEG result |224| on 
the /i — > e7 branching ratio, Br(fi — > ey) < 2.4 x 10 -12 at 95% C.L. and to other similar processes 
like r — > (e or fi)j. One expects that lepton flavour-violating processes may have a large discriminating 
power in assessing the relative merits of the different models proposed for neutrino mixing. In fact, one 
must pay attention that the large corrective terms introduced to shift 6*13 from the TB value could appear 
in the non-diagonal elements of the charged lepton and s-lepton mass matrices (in a basis where all kinetic 
terms are canonical) and could induce a too large // — > ej branching ratio ]225f [230|. This problem has 
been discussed in detail in ref. (58) within the simple CMSSM framework (Constrained MSSM). While 
this GUT-constrained version of supersymmetry is rather marginal after the results of the LHC searches, 
more so if the Higgs mass is confirmed to lay around to# = 125 GeV, still we think that it can be used for 
our purposes in the present context. 

The results derived in ref. [58 1 (see|5ja)) show that the typical A4 models are well suited to satisfy 
the MEG experimental bound, as the non diagonal charged (s-)lepton matrix elements needed to best 
approximate the mixing angles are particularly small in these models. In fact, the size of the dangerous off 
diagonal terms is driven by the value of the parameter £ whose optimal value in the procedure of sect. 4 was 
found to be £ ~ 0.07. This value corresponds to a modest score in terms of success rate for reproducing 
the mixing angles in a scanning of the parameter space ( ~ 18% ), but is sufficiently small to maintain 
the off diagonal charged (s-)lepton mass matrix elements within affordable limits given the bounds on 
lepton flavour violating processes. A comparable score is also achieved by the models of the Lin type (see 
see |5jb)), because the (on the average) larger correction to the mixing angle 613 actually arises from the 
neutrino sector in these models, while the corrections from the charged lepton sector are naturally kept at 
a smaller level. 
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Mi/2 (GeV) M m (GeV) 

(a) Typical A 4 : tan = 2 and m = 200GeV (b) Special j4 4 : tan /3 = 2 and m = 200GeV 

Fig. 5 Scatter plots of BR(/i — > ej) as a function of M1/2, for tan j3 — 2, and mo = 200 GeV |58| . The parameters 
£ (for typical A4) and £' (for special A4) are chosen in order to maximize the success rate of each model (£ = 0.075, 
|£'| = 0.183). For Blue (Red) points the lightest supersymmetric mass (LSP) is the lightest neutralino (stau). The 
percentage in each plot refers to the number of Blue points that satisfy the MEG bound over the total number of 
points. The horizontal line shows the current MEG bound. For larger values of tan /3 the success rate decreases, while 
increases for larger values of mo, being more difficult for the stau to be the LSP. 

In conclusion, when the fit to the mixing angles and the bounds on lepton flavour violating processes 
are combined, the typical A4 models are rather weak on the mixing angles but, as discussed in detail in 
Ref. [58 1, are better suited to cope with the bounds on lepton flavour violating processes. The special A 4 
models of the Lin type offer the best overall performance to the data as they are far better on the mixing 
angles and comparable, although at a lower level, on lepton flavour violating processes. As for the regions 
of the CMSSM parameter space that are indicated by our analysis the preference is for small tan/3 and 
large SUSY masses (at least one out of mo and mi/2 must be above 1 TeV). As a consequence it appears 
impossible, in these models, at least within the CMSSM framework, to satisfy the MEG bound and, at the 
same time, to reproduce the muon g — 2 discrepancy. 



9 Conclusion 

The recent rather precise measurements of #13 make our present knowledge of the neutrino mixing matrix, 
except for the CP violating phases, sufficiently complete to considerably restrict the class of models that 
can reproduce the data. In spite of this process the range of possibilities remains unfortunately quite wide. 
On the one extreme, the rather large value measured for (9 13 , close to the old CHOOZ bound, has validated 
the prediction of models based on anarchy [231 232], i.e. no symmetry in the leptonic sector, only chance, 



so that this possibility remains valid, as discussed, for example, in ref. 1 233 1. Anarchy can be formulated 
in a SU (5) U (1) context by taking different Froggatt-Nielsen J98) charges only for the SU(5) tenplets 
(for example 10 ~ (3, 2, 0), where 3 is the charge of the first generation, 2 of the second, zero of the third) 
while no charge differences appear in the 5 (5 ~ (0, 0, 0)). Anarchy can be mitigated by assuming that 
it only holds in the 2-3 sector with the advantage that the first generation masses and the angle 6*13 are 
naturally small (see also the recent revisiting in ref. (234)). In models with See-Saw, one can also play 
with the charges for the right-handed SU(5) singlet neutrinos. If, for example, one takes 1 ~ (1, —1, 0), 
together with 5 ~ (2, 0, 0), it is then possible to get a normal hierarchy model with #13 small and also with 
r = Aml olar / Am^ tm naturally small (see, for example, Ref. [235 1). In summary anarchy and its variants, 
all based on chance, offer a rather economical class of models that are among those encouraged by the new 
813 result. On the other extreme, stimulated by the fact that the data suggest some special mixing patterns 
as good first approximations, in particular TB mixing, models based on discrete flavour symmetries, like 
A4 or 54, have been proposed and widely studied. 
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In A4 models, the A4 symmetry is broken down to two different subgroups in the charged lepton sector 
and in the neutrino sector, and the mixing matrix arises from the mismatch between the two different 
residual symmetries. The breaking can be realized in a natural way through the specific vacuum alignments 
of a set of scalar flavons. There are many variants of models where TB mixing is indeed derived at leading 
order (in particular with or without see-saw) with different detailed predictions for the spectrum of neutrino 
masses and for deviations from the TB values of the mixing angles. The starting LO approximation is 
completely fixed (no chance), but the NLO corrections introduce a number of undetermined parameters. In 
general at NLO the different mixing angles receive corrections of the same order of magnitude, which are 
constrained to be small due to the experimental results which are close to the TB values. Indeed the small 
experimental error on (9 12 , with a central value that is close to the value predicted by TB mixing, suggests 
that the NLO corrections should be of order of few percent, at most. The recent data on 13 and the MEG 
new upper bound on the LFV process /1 —> ej impose a reappraisal of these models [58]. In particular, the 
relatively large value of 613 introduces a marked departure from the TB limit, while the values of 6*12 and 
823 are very close to it. The challenge is to produce in a natural way a relatively large correction to #13 
without affecting too much the other mixing angles. But one must pay attention that these larger corrective 
terms introduced to shift (9 13 from the TB value could appear in the non-diagonal elements of the charged 
lepton (and s-lepton) mass matrix and could induce a too large \i — > e-f branching ratio. 

As a result of a detailed analysis p8[ we find that, for reproducing the mixing angles, the Lin type A4 
models have the best performance, as expected, but the typical A4 models can also accommodate the data 
with a reasonable probability. As for lepton flavor violating processes, the problem has been studied by 
adopting the simple CMSSM framework. While this over constrained version of supersymmetry is rather 
marginal after the results of the LHC searches, more so if the Higgs mass really is around mjj = 125 
GeV, still we think it can be used here for indicative purposes. The typical A4 models turn out to be the 
best suited to satisfy the MEG experimental bound, as the non diagonal charged lepton matrix elements 
needed to reproduce the mixing angles are rather small. A slightly worse score, but still rather good, is 
achieved by the models of the Lin type, where the main corrections to the mixing angles arise from the 
neutrino sector. When the fit to the mixing angles and the bounds on LFV processes are combined, the A 4 
models emerge well from our analysis and in particular those of the Lin type are remarkably successful 
in the lepton sector. As for the regions of the CMSSM parameter space that are indicated by our analysis 
the preference is for small tan/3 and large SUSY masses (at least one out of mo and mi/2 must be above 
1 TeV). As a consequence it appears impossible, at least within the CMSSM model, to satisfy the MEG 
bound and, at the same time, to reproduce the muon g — 2 discrepancy. 

It is remarkable that neutrino and quark mixings have such a different qualitative pattern. An obvious 
question is whether some additional indication for discrete flavour groups can be obtained by considering 
the extension of the models to the quark sector, perhaps in a Grand Unified context. The answer appears 
to be that, while the quark masses and mixings can indeed be reproduced in models where TB (or GR or 
BM) mixing is realized in the leptonic sector through the action of discrete groups, there are no specific 
additional hints in favor of discrete groups that come from the quark sector fl4) . 

Finally, one could have imagined that neutrinos would bring a decisive boost towards the formulation 
of a comprehensive understanding of fermion masses and mixings. In reality it is frustrating that no real 
illumination was sparked on the problem of flavour. We can reproduce the observations in many different 
ways, in a wide range of models that goes from anarchy to discrete flavour symmetries, but we have not 
yet been able to single out a unique and convincing baseline for the understanding of fermion masses and 
mixings. In spite of many interesting ideas and the formulation of many elegant models the mysteries of 
the flavour structure of the three generations of fermions have not been much unveiled. 
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